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The investigation of the properties of a Higgs boson, especially a test of the predicted 
linear dependence of the branching ratios on the mass of the final state, is currently 
one of the most compelling arguments for building a linear collider. We demonstrate 
that the large Higgs boson production cross section at a 3 TeV CLIC machine allows 
for a precision measurement of the Higgs branching ratios. The cross section times 
branching ratio of the decays H — > 66, H — >^ cc and H — > can be measured with a 
statistical uncertainty of 0.22%, 3.2% and 15%, respectively. 

The Higgs mechanism of clcctrowcak symmetry breaking predicts the existence of a 
fundamental spin-0 particle that has not yet been discovered. The LHC will most probably 
answer the question about its existence by 2013. The Standard Model predicts a linear 
dependence of the Higgs branching ratios on the mass of the final state, but non-standard 
couplings could alter this relation. The detailed exploration of the Higgs sector will be 
instrumental to our understanding of the fundamental interactions. The LHC will be of 
limited utility in this exploration. 

The compact linear collider (CLIC) is a proposed c"'"e~ collider with a maximum center- 
of-momentum energy ^/s = 3 TeV, based on a two-beam acceleration scheme [1]. We present 
here the analysis of the branching ratios H — > 66, H — > cc [8] and H — )■ [j.^ [4] at such a 
machine. The studies are based on fully simulated samples in the CLIC_SiD [5] detector 
concept and take into account the main beam-related background. 

1 The CLIC_SiD Detector Model 

The detector that is used in the full simulation of samples is based on the SiD concept [2] 
and has been adapted [5] to the specific detector requirements at CLIC. It is designed for 
particle flow calorimetry using highly granular calorimeters. 

A superconducting solenoid with an inner radius of 2.9 m provides a central magnetic 
field of 5T. The calorimeters are placed inside the coil and consist of a 30 layer tungsten- 
silicon electromagnetic calorimeter with 3.5 x 3.5 mm^ segmentation, followed by a tungsten- 
scintillator liadronic calorimeter with 75 layers in the barrel region and a steel-scintillator 
hadronic calorimeter with 60 layers in the endcaps. The read-out cell size in the hadronic 
calorimeters is 30 x 30 mm^. The iron return yoke outside of the coil is instrumented with 
9 double RPC layers with 30 x 30 mm^ read-out cells for muon identification. 

The silicon-only tracking system consists of 5 20 x 20 ^m^ pixel layers followed by 5 strip 
layers with a pitch of 25 |im, a read-out pitch of 50 |a,m and a length of 92 mm in the barrel 
region. The tracking system in the endcap consists of 5 strip disks with similar pitch and a 
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Table 1: List of processes considered for this analysis with their respective cross section a 
and the number of simulated events A/events- The cross section takes into account the CLIC 
luminosity spectrum. Cross sections marked with * include a cut on the invariant mass of 
the muon pair to lie between 100 and 140 GeV. 

stereo angle of 12°, complemented by 7 pixelated disks in the vertex and far- forward region 
at lower radii with pixel sizes of of 20 x 20 (xm^. 

2 Softweire and Data Samples 

The physical processes arc produced with the Whizard [7,9] event generator, with fragmen- 
tation and hadronization by PYTHIA [10]. The events are simulated in the CLIC_SID 
detector model using SLIC, a thin wrapper around GEANT4 [3]. The event reconstruction 
is handled by the Icsim and slicPandora packages, the LCFI package is used for flavor tag- 
ging. The assumed luminosity of the analyses is 2 ab~^, corresponding to about 4 years of 
data taking at nominal conditions. 

The predominant production channel of Higgs bosons at a 3 TeV linear collider is the 
WW fusion channel c+e^ — > hvv. The main backgrounds arc from events with two Z 
bosons, where one Z decays to the signal final states, and the other decays to neutrinos, 
or decays to a fermion pair that does not enter the fiducial volume of the detector. The 
beam configuration at a 3 TeV CLIC machine produces 3.2 yy hadrons events per bunch 
crossing on average. With a spacing of 0.5 ns between bunches, these necessarily pile up 
in the subdetectors, which have integration times of 10 ns, except for the barrel hadronic 
calorimeter, which has an integration time of 100 ns. To take into account the effect of this 
background on the measurement, a fully simulated sample of events from yy — > hadrons 
corresponding to 60 bunch crossings is mixed with each physics event for the analysis of 
the Higgs decaying to b and c quarks. In the H — >■ (ip, analysis, only the signal sample was 
mixed with events from yy — > hadrons background. Table 1 lists the physics processes that 
were taken into account in the analyses, together with their cross section and the number 
of simulated events. 
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3 Measurement of H — > 66 and H — > cc 



3.1 Event Selection 

The main background of the measurement of the decays H — )• 66 and H cc is from two- 
jet processes e+e" qqw, due to their large cross section, and from processes with two 
measured jets and additional particles that escape detection. 

The Fast Jet package is used to cluster the events into two jets. The LCFI flavor tagging 
package finds secondary vertices in each jet and uses them in a neural network to distinguish 
b-, C-, and light quark jets. The invariant mass of the jet pair is the major discriminant 
between decays of Higgs and of Z bosons. It is used in a second neural network, together 
with the output of the flavor-tagging network and the following variables: 

• the maximum of the absolute values of jet pseudorapidities 

• the sum of the remaining LCFI jet flavor tag values, i.e., c(udsb), c(b)-tags and b(uds)- 
tag 

• Rr,(t>, the distance of jets in the r/ — plane 

• the sum of jet energies 

• the total number of leptons in an event 

• the total number of photons in an event 

• the acoplanarity of the jets 

The neural network selection efiiciency S/{S + B) versus the statistical uncertainty 
\/S + B/S on the measurement of the number of signal events S and background events B 
is shown in Figure 1 for the two neural networks that were trained on H ^ 66 and H ^ cc 
as signal, respectively. The optimal selection is at the local minimum of the curve, at a 
selection efiiciency of 55% for H — >■ 66 with a sample purity of 65%, corresponding to a 
statistical uncertainty of 0.22%. The optimal selection for H — > cc has an efficiency of 15%, 
corresponding to a sample purity of 24% and a measurement uncertainty of 3.2%. These 
values reflect the fact that b-jets can be distinguished from c-jets with high purity, while in- 
completely reconstructed b-jets make up a large fraction of the background to c-jet selection, 
making the analysis more challenging. Using the output of the reconstruction algorithms in 
neural networks leads to the minimal statistical uncertainty on the measurement at the cost 
of an increased dependence on systematic effects. We assume that with sufficient expericince 
at the running machine, the systematic variations are well enough understood so that the 
systematic uncertainties are comparable to the statistical uncertainties. 

4 Measurement of H — > 

The measurement of the rare decay H — >• ma, requires high luminosity operation and sets 
stringent limits on the momentum resolution of the tracking detectors. The branching ratio 
of the decay of a Standard Model Higgs boson to a pair of muons is important as the lower 
end of the accessible decays and defines the endpoint of the test of the predicted linear 
dependence of the branching ratios to the mass of the final state particles. 
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Selection efficiency Selection efficiency 

(a) H -> bS (b) H -> cc 

Figure 1: Statistical uncertainty of the measurement of cross section times branching ratio 
versus selection efficiency of the neural network. Left: The neural network was trained 
to identify H — )• 66 decays from di-jet backgrounds including H — )■ cc. Right: The neural 
network was trained on H —>■ cc as signal and di-jets including H — > 66 as background. 
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The muon reconstruction efhciency in the 
signal sample is shown in Figure 2. The 
small deterioration of the muon reconstruc- 
tion efficiency due to beam-induced back- 
ground from YY hadrons is evident. The 
average muon reconstruction efficiency for 
polar angles greater than 10° is 98.4% with 
this background compared to 99.6% with- 
out. The total reconstruction efficiency 
of the signal sample, requiring two recon- 
structed muons with an invariant mass be- 
tween 105 GeV and 135 GeV is 72% in 
the presence of background. The events 
are selected by requiring two reconstructed 
muons, each with a transverse momentum 
of at least 5 GeV. In case there are more 
than two muons reconstructed, the two 
most energetic ones are used. In this note, 
the most energetic muon is referred to as (ii 
and the second most energetic muon is re- 
ferred to as ii2- In addition, the invariant mass of the two muons M{\iii) is required to lie 
between 105 GeV and 135 GeV. 
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Figure 2: Muon reconstruction efhciency for 
the signal sample with and without yy — >■ 
hadrons pile-up 
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The event selection is done using the boosted decision tree classifier implemented in 



TMVA [6]. The 



and x+t vv samples are not used in the training of the BDT, 



but are effectively removed by the classifier nevertheless. 

The variables used for the event selection by the BDT are: 

• The visible energy excluding the two reconstructed muons -Evis- 

• The scalar sum of the transverse momenta of the two muons pt(^i) +pt(|J-2)- 



The helicity angle cos 6** (fin) 



P (Hl)-P(n^) 



where is the momentum in the rest 



|p'(m)|-|p(tiH)l ■ 

frame of the di-muon system. Since the two muons are back-to-back in the rest frame 
of the di-muon system there is no additional information to be gained from calculating 
a similar angle for \x.2. 

• The relativistic velocity of the di-muon system P(|j,(i), where P = ^. 

• The transverse momentum of the di-muon system pt(|J-M-)- 

• The polar angle of the di-muon system 9{\m). 

The most powerful variable is the visible energy whenever there is an electron within the 
detector acceptance. Otherwise the background can be rejected by the transverse momentum 
of the di-muon system or the sum of the two individual transverse momenta. Figure 3 
clearly shows the Higgs peak in the invariant mass distribution after the event selection. 
The background from e+e" — )■ e"'"e~(i|j. 
events, is effectively reduced by forward 
electron tagging. While the forward 
calorimeters were not part of the full de- 
tector simulation, assuming a tagging effi- 
ciency of 95% down to an angle of 40 mrad 
for electrons of several hundred GeV to over 
one TeV is a conservative estimate, even in 
the presence of yy — > hadrons background. 
It is found that Bhabha events prevent fur- 
ther rejection of this background at lower 
angles. The results quoted are based on a 
ad-hoc rejection of 95% of the electrons in 
the Luminosity Calorimeter. 

5.1 Invariant Mass Fit 
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Figure 3: Maximum Likelihood fit of the Higgs 
mass in the data sample after selection cuts 



The distribution of the invariant mass in the 
H — >■ |in sample has a tail towards lower 
masses because of final state radiation. The 
shape can be described best by two half 

Gaussian distributions with an exponential tail. Together with the mean value this results 
in five free parameters in the fitted function, which can be written as 



g2.£+ai(.-™o)- , a; < mo 



, X > mo 



5 



LCWSll 



where mo is the mean of both Gaussian distributions, and an are the widths, and 
and aR arc the tail parameters of the left and the right Gaussian distribution, respectively; 
n is a normalization parameter. The background is well described by an exponential param- 
eterization, obtained from a background-only sample. 

The number of signal events is obtained from a maximum likelihood fit to the sample 
containing signal plus background after the event selection. 



5.2 Study of the Momentum Resolution 



The ability to measure the decay H — > (X|j, 
depends crucially on the momentum reso- 
lution of the tracking detectors. In a fast 
simulation study, different values for the 
momentum resolution were assumed, before 
the sample was fit as described in the pre- 
vious section. The results are shown in Ta- 
ble 2. It is found that an average resolution 
of 5x 10~^ GeV~^ or better is required in or- 
der for the momentum resolution not to be 
the dominant uncertainty contribution in a 
2ab~^ measurement of the decay H 
The average muon momentum resolution of 
the fully simulated sample is 4x 10~^ GeV~^ 
corresponding to a statistical uncertainty of 
23% without forward electron tagging. 

5.3 Results 
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Table 2: 



h — > |J,+|J, 



Summary of the results for the 

branching ratio measurement us- 
ing fast simulation samples with different mo- 
mentum resolutions cr (Apt)/Pt assuming an 
integrated luminosity of 2ab-^. Given is 
the corresponding invariant mass resolution 
(T(AM(|i|a,)) and the resulting statistical un- 
certainty of the (7hv Vo X BR-h^n+|i- measure- 
ment. The results do not include reduction of 



the e+e 



H|x background using ad-hoc 



electron tagging. 

We have demonstrated the feasibility of 
measuring the branching ratios of a 120 
GeV Higgs boson at a 3 TeV CLIC with 
high precision. For the measurement of 

Higgs decays to quarks, 0.22% and 3. 2%. statistical uncertainty can be achieved for the 
decays H — > 66 and H — >• cc, respectively. This includes the effect of background from 
Y/ hadrons on the flavor tagging. 

For the rare decay H the cross section times branching ratio can be measured to 



a precision of 15% if the background from e+e 



can be reduced using tagging of 



electrons down to an angle of 40 mrad with an efficiency of 95%, and the average momentum 
resolution is not worse than 5 x 10"®. The effect of background from yy hadrons has been 
taken into account. From experience of the LEP experiments one can assume that the 
systematic uncertainties related to detector effects are of the order of 1% or less. For the 
measurement of <JTfi^\i+v.- at LEP the systematic uncertainty was between 0.1 and 0.4%, 
depending on the experiment. In summary, one can expect that the systematic uncertainty 
of this analysis will be negligible compared to the statistical uncertainty. 

The expected uncertainty of the peak luminosity is currently being studied but is esti- 
mated to be around 1% or less. 
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